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AilTRACT 


A  dl/fuslon- deposition  nxiel  Is  developed  for  use  In  estimating 
dosage  levels  duo  to  in-r flight  lease  of  fission  fragments  from  a 
nuclesu'-powere^  >  aircraft.  T>*«  model  is  based  on  the  work  of  Sutton 
for  dUfusioc  <uid  the  work  of  Coamberlaln  for  deposition.  The  model 
eonsiders  an  flevnted  kistar-l^  jous  point  source  and  an  elevated 
Inatantaneoas  line  source  oriented  at  an  arbitraiy  angle  to  the  mean 
wind  (Urectioii, 

Suggested  values  of  the  various  dsposltiony  rain- out.  and  diffusion 
parameters  to  be  used  with  the  model  are  presented  along  with  a  quail- 
tative  dlscusalon  of  the  uncertatntisa  of  the  model  and  the  suggested 

parameters. 
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/:  DIFFUSION-DEPOSITION  MODEL 
FOR  IN-FLIGHT  RELEASE  OF  FISSION  FRAGMENTS 

1.  INTRODUCTION 

The  purposes  of  this  study  have  been  (1)  to  develop  a  diffusion- 
deposition  model  which,  when  combined  with  allowable  fission- product 
dosages,  will  permit  estimation  of  the  allowable  fission- product  release 
rates  during  flight  operations  of  nuclear  aircraft,  and  (2)  to  suggest 
values  of  the  parameters  to  be  used  in  the  model  for  a  variety  of 
meteorological  conditions,  flight  altitudes,  and  flight  patterns. 

As  a  first  approximation  and  with  no  slgniUicant  loss  of  generality, 
the  source  of  contamination  has  been  treated  either  as  a  line  source  of 
finite  length  or  as  a  point  source  at  each  of  several  fixed  heights  above 
sea  level.  Values  of  the  parameters  have  been  selected  for  flight  alti¬ 
tudes  of  nuclear  aircraft  ranging  from  500  to  35, 000  feet  above  sea 
level.  Flights  in  the  stratosphere  have  not  been  considered. 

It  has  been  assumed  that  the  fission  products  will  be  released  as 
gases  or  as  fine  particulates.  The  particles  are  expected  to  exhibit  a 
log  normal  distribution  with  a  geometric  mean  diameter  of  0.03^  and  a 
geometric  standard  deviation  of  two.  Also  assumed  is  a  constant  release 
rata  of  each  fission  product  during  periods  of  operation  on  nuclear  fuel. 

Consideration  has  been  given  to  the  problem  of  exposure  not  only 
of  the  fixed  receptor  at  ground  level  but  of  the  moving  receptor  at  air¬ 
craft  flight  levels. 

2.  THE  DIFFUSION  MODEL 

2. 1  Definition  of  Problem;  Modeling  Approximations 

The  diffusion  model  is  derived  to  apply  to  the  following  specific 

problems  associated  with  the  normal  operating  release  of  radioactive 

materials  by  a  nuclear- powered  aircraft; 

a.  A  source  geometry  of  either  a  finite  straight  line  or 
a  circle  laid  out  at  some  given  height  between  500  unci  35,000 
ft  above  the  ground  surface 

(Authors'  manuscript  approved  13  April  lOGO.) 


1).  An  orientation  of  the  line  source  at  an  arbitrary 
aJ'Jtle  to  the  mean  wind  ciire;tlon  at  lll^ht  altitude 

c.  An  estimate  of  the  i.'ollution  levels  to  which  a 
person  on  the  ground  may  be  exposed 

d.  An  estimate  of  the  pollution  levels  at  an  avbj'.rary 
heigiit  in  the  atmosphere  to  which  a  person  flying  iii  a 
private  or  commercial  aircraft  may  be  exposed 

The  finite  time  the  aircraft  will  take  tn  lay  out  the  source  has  been 
ignored  In  development  of  the  diffusion  models.  That  Is,  It  has  been  as¬ 
sumed  that  the  source  may  be  considered  Instantaneous  in  comparison 
with  the  diffusion  times  that  characterize  the  problem  (several  hours  to 
days).  This  reasoning  also  leads  to  the  assumption  that  the  source  laid 
out  In  a  circle  may  be  treated  as  an  instantaneous  point  source  for  ex¬ 
posure  distances  sufficiently  far  from  the  circle.  The  error  introduced 
by  these  approximations  Is  effectively  nil  for  estimating  contamination 
levels  at  the  ground  and  at  aircraft  flight  levels. 

The  must  serious  assumption  made  In  the  models  is  that  the  mean 
wind  Is  constant  in  direction  and  speed  throughout  the  portion  of  the  at¬ 
mosphere  under  consideration.  In  general  this  is  an  unrealistic 
assuHiptlon  for  application  to  "typical"  atmospheric  conditions.  Never¬ 
theless  It  Is  felt  that  the  uncertainties  arising  In  the  selscticn  of  values 
for  the  various  modeling  parameters  arc  large  enough  that  trying  to 
account  rigorously  for  a  nonuniform  mean  wind  field  introduces  an 
unnecessary  refinement.  A  qualitative  discussion  of  the  error  this 
assumption  Introduces  is  given  In  Section  3  of  tlUs  report. 

The  basic  diffusion  model  chosen  as  a  starting  point  for  this 
12 

problem  Is  Sutton's  model  for  concentrations  due  to  an  instantaneous 
poliit  source.  The  final  model  developed  for  the  ground  contamination 
problem  enables  one  to  malce  rough  estimates  of  the  following  quantities 
for  either  the  point  source  or  the  line  source: 

a.  The  exposure  or  tlme-lntcgrated-concentratlon  with 
or  without  deposition  on  the  ground  or  scavenging  ty  pre¬ 
cipitation  processes 


2 


b.  The  total  deposition  of  material  at  a  point  on  the 
ground  due  to  natural  pickup  of  the  material  l)y  vegetation, 
objects,  etc. 

c.  The  total  deposition  of  material  at  a  poLnt  on  the 
ground  due  to  scavenging  by  precipitation  processes 

The  model  for  estimating  exposures  at  flight  altitudes  ignores  the 
effect  of  deposition  or  scavenging  on  the  air  conceidrations.  Only  two 
special  cases  of  this  problem  are  treated:  an  airplane  flying  through 
the  cloud  fl)  parallel  to  the  flight  line  of  the  nuclear-powered  aircraft 
at  the  same  speed,  and  (2)  parallel  to  the  mean  wind  direction. 

Finally,  In  both  models  li  is  assumed  that  there  is  no  difference 
between  the  power  Indices  on  distance  for  diffusion  in  the  crosswind  and 
vertical  directions:  tiiat  Is,  It  is  aisumed  that  any  anisotropy  existing 
In  the  atmosphere  on  the  scale  of  this  problem  can  be  accounted  for  by 
simply  assigning  different  values  to  the  so-called  diffusion  coefficients, 
Cy  and  Cj  , 

2, 2  Model  for  Estimating  Contamination  at  Ground  Level 

The  equation  for  the  concentration,  x  •  ^  on.  the  grouno, 

(x,  y,  0),  due  to  an  Instajntaneous  point  source  located  at  x  a  /cos  0, 
y  a  .I'sln  Q,  z  a  h  (see  Fig.  1)  Is  given  by 


X  (*.y.o.O  ^ 


2Q* 
X  y  *  • 


exp  (  - 


1 

(x-^cosO  -ut)^  (y-^'sln  ?)^ 

1 

C  ^  C  ^  ^ 

L  X  '"y  '^zJ 

(1) 


where  Q*  Is  the  source  strength  (total  amount  of  material);  ,  C 

and  C-  are  the  diffusion  coefficients  with  the  dimensions  of  length  to 
n  * 

the-y  power;  n  is  the  power  index  (dlmen  lonless)  mentioned  prev¬ 
iously;  V  has  the  dimensions  of  quantity  of  material  per  unit  volume; 


3 


1«  Schematic  diagram  of  ceometrv  at 

dir^^fnr,^?  parallel  to  the  mean  wind 

curectlon.  Length  of  source  »  2L;  angle  between 
line  source  and  x-axls  »  9;  height  ouSw 
above  ground  plane  (z  =f  o) h. 

u  is  the  mean  wind  speed;  t  la  time  after  introduction  of  the  source- 

nd  X.  y  and  z  are  the  coordinates  of  a  Cartesian  system  with  the  ’ 
x-a;as  oriented  parallel  to  the  mean  wind  direction  (Fig.  i).  The 
origin  of  the  coordinate  system  is  the  point  defined  by  i  «  o  (x  =:  y  o), 

The  total  exposure,  E  .  is  obtained  from  Eq.  (1)  by  Integratlnir 
over  time  (zero  to  Infinity)  assuming  that  the  spread  of  the  cloud  Is 

Cloud  fromT'l  .h  T 

Cloud  from  its  Initial  posUlon,  The  result  Is 


1 


,  t'<,  y.  'V' 


y..  C  ,  1  f  x  c;j  .  ;  -  ' 


(y-^sln-)"  h 


exp  s  - - ; - — -  - - - r  — ^5- 

{X  - 'COS  I  c  2 


Th«i  aruiogcus  quantity,  £,  for  an  Instantaneous  line  source  is 
obtained  by  Integrating  Eq.  (2)  over /;  from  -  ^  to  +L  for  a  line 
source  of  length  2L: 


2Q  /,(x)  I 

E(x,y.  0)  a  -  / - 

rtCyCjil  •/l  (x cosO)'^-’' 


(X  - /cos  ;)2-n 


(y  -  /  sin  ""I 

- ^  (3) 


where  b(x)' 


for  X  ''  L  cos  '•  ;  ^  00’ 


L  for  X  >  L  cos  3 


L  for  •  -  00®  :  X  >  o 


The  source  strength,  Q  ,  must  now  be  given  in  terms  of  amount 
of  material  per  unit  length  to  maJce  the  equation  dimensionally  con¬ 
sistent. 

If  there  were  no  deposition  or  precipitation  scavenging,  Eq.  (3) 
would  give  estimates  of  the  exposure  near  the  grou;ul  upon  selection  of 
suitable  values  for  the  various  parameters.  It  Is  felt,  however,  that 
these  latter  processes  may  be  Important,  particularly  for  thnse  re¬ 
leases  relatively  near  the  ground  surface.  An  atten^pt  to  account  for 
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these  effects  Is  made  following  the  work  of  Chamirer'aln. 

The  total  deposition,  2.  ,  at  a  point  on  the  grour.d  Is  given  by 


'Q 

-(X.  y.  o>  =  Vv  E(x,  y,  o)  =  -  —  (• 

where  has  the  dimensions  of  quantity  of  material  per  unit  area; 

V\  Is  the  deposition  rate.  The  far  right-hand  side  of  Eq,  (5)  indicates 
S 

the  fact  that  deposition  acts  to  deplete  the  effective  source  strength  of 
the  cloud.  Substituting  Eq.  (3)  into  Eq.  (5),  solving  for  Q,  and  sub¬ 
stituting  the  result  back  into  Eq.  (3)  gives 

E(x,y,o)  =*  -  jexp  -  /  / - r—- 

^C  C  u[  '-LcoaO’'-L  tC  C.iT (x -c  cos 
y  z  y  z 


exp 


1 


(X-;fCO8  0)2**'  Cy2 


/(y-.  sin 9)2  h' 


'V 


d>-dx  1 


(0) 


^  / 


..i(x) 


1 


-L  (x -icos  9)2-*' 


exp  - 


(x cos  r)2-*i  \  Cy2 


(ywsin")"  h^ 


d  ! 

i 

) 


where  is  the  source  strength  at  x  a  o  (or  at  the  initial  time).  The 
deposition,  ^  :  is  then  given  by  multiplying  Eq.  (6)  by  Vg  ,  the 
deposition  rate. 

The  scavenging  of  the  particulate  or  contaminant  cloud  by  pre¬ 
cipitation  processes  is  treated  here  as  a  two- stage  process.  Decause 
of  the  small  particle  sizes  expected  in  tho  contaminate  cloud,  one  must 
first  estimate  the  fraction  of  exposure  that  becomes  contained  in  the 
cloud  droplets.  The  second  step  is  to  estimate  the  efficiency  of  pre¬ 
cipitation  (rain)  in  removing  the  cloud  droplets.  The  actual  precipita¬ 
tion  depletes  the  effective  source  strength  of  the  particulate  cloud 
whereas  the  acquisition  of  particles  by  cloud  droplets  has  no  effect  on 
the  exposures.  Once  precipitation  has  started,  the  exposure  at  the 
ground  (following  Chamberlain^)  is  given  by 
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E(x.y,c)  = 


1 


where  A  Is  the  elimination  constant  (units  of  time'  that  characterizes 
the  depletion  of  cloud  droplets  by  precipitation;  Is  the  downwind 
distance  where  precipitation  scavenging  first  begins;  that  is, 


I"  i  «  < 

,  A  ;  X  k  Xjj  .  (8) 


la  the  time  constant  describing  the  cloud  droplet  efficiency  in  acquiring 
particles,  and  x_  is  the  downwind  distance  where  cloud  droplets  first 

A 

appear.  The  total  amount  of  material  per  unit  area  deposited  on  the 
ground  is  obtained  by  integrating  >'\.  [l  >  exp  ^  ^ 

A 

over  height  under  the  assumption;  that  the  removal  processes  are 
uniform  with  height  and  that  the  water-droplet  cloud  completely  em- 
I)race3  the  particulate  cloud.  Dry  deposition  and  reflection  at  the 
^surface  when  precipitation  occurs  are  both  Ignored.  Thus 


P(x,y,o)  =»  1  -  exp  -  :  (X  -  x  )/u  ■,  ,  ‘  E{x,y,  2)  dz 

i  '0 

I  J' Qo '-'xp  [-  (x-Xb)/ilj 

=  a-exp  -v(x-x^)/i:^|--  . 


i(jt)  1 

.-L  (x-xCoa9)^ 


'  (y  -  ^  sLnOr  j 
[  C  2  (X  .  /cos  oT^j 


If  a  poLnt  source  of  streiitith  Instead  of  a  finite  line  source  la 
considered,  the  above  modelii^  equations  are  simplified  by  setting 
identically  equal  to  zero  v;hich  eliminates  the  inlegrati«3Rs  with  respect 
to  .c  .  Thus,  to  summarize,  tne  solutions  for  the  point  source 
problem  are  as  fallows: 

a.  Exposure,  with  dry  deposition. 


Ep(x.y.o) 


•  0 

cy 

1 

i  1 

/  y2 

h^ 

A. - 4 

:  C,  tix2-n 

1  ■*  On 

(  x2-'' 

C  2  / 

y  * 

y 

4 

r 

1  X 

:  '1 

2Vg 

*w 

i  1 

/  1 
y' 

exp  •  -  /  — 


Jo  TiC  C  U  X"'**' 

y  ^ 


,2-n  'p  2  r  ^  / 


b.  Deposition 
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c.  Exposure,  .vith  precipitation  scavenging, 
2Qj  exp  C-/,(x  -  xjj)/a  j 


E„(x,  y,  o)  =  : 


TtCyC^u  x2-n 


exp 


(13) 


d.  Deposition  due  to  precipitation  scavenging, 


Pp(x.y.o)  =» 


r  i!  A Qo*  C-A(*  -  Xb)/^  3 


exp 


o  2  _2-  n 

S  * 


(14) 


2. 3  Model  for  Estimating  Contamination  for  Privately  Flown  Aircraft 

Let  us  Imagine  a  privately-owned  or  commercial  plane  flying  a 
course  parallel  to  the  source  line  laid  out  by  the  nuclear- powered  air¬ 
craft,  To  simplify  the  problem  mathematically,  it  is  assumed  that  the 
planes  have  equal  flight  speeds  and  that  the  end  effects  due  to  Uie  finite 
length  of  the  source  can  be  neglected.  The  source  can  now  be  treated 
as  a  fixed  continuous  point  source  axul  the  exposed  plane  as  a  fixed  point 
relative  to  the  scurce.  Surface  deposition  does  not  enter  this  problem 
and  precipitation  scavenging  is  neglected.  The  conrentration  at  the 
exposed  plane  is  then  given  by 


X(x,y.z) 


!i  ! 
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wiiere  the  Cartesian  coordinate  systems  has  Its  origin  at  the  exhaust 
point  of  the  nuclear -powered  aircraft;  Q  Is  the  amount  of  material 
released  per  unit  time. 

The  exposure  for  the  private  plane  Is  given,  within  the  l‘'nit3  of 
the  preceding  assumptions,  by 


E(x,y.*) 


-C  C  x^-'^ 
/  z 


—  exp 


(16) 


where  Q  Is  the  total  amount  of  material  released  In  flying  a  line  ut 
length  2L;  that  U 

Q  5*  /  Q  dt 


where  VL  Is  the  airspeed  of  the  airplane.  This  case  is  included  be* 

*X 

cause  U  serves  to  define  a  region  in  the  vicinity  of  the  nuclear -powered 
aircraft  from  which  other  air  traffic  must  be  restricted.  It  is  the  most 
serious  conceivable  air  contamination  problem  arising  from  normal 
operation  of  the  nuclear -powered  aircraft. 

The  case  of  a  privately-owned  plane  flying  parallel  to  the  mean 
wind  direction  is  treated  here  for  an  infinite  line  source  oriented  per¬ 
pendicular  to  the  mean  wind.  This  is  not  a  particularly  serious  simpli¬ 
fication  if  the  exposed  plane  flies  along  the  x-axls  of  the  finite  lino 
source;  in  effect,  the  model  neglects  the  end  effects  of  the  real  finite- 
length  line  source. 

The  concentration  at  a  point  (x,  z)  relative  to  the  initial  release 
line  is  given  by 


;<(x,z,t)  =1 


Qo 

<■ 

1 

*(x-ut)* 

■  exp 

:tCyC^(ut)2-« 

{ut)2-n 

p  2  r*  ^ 

^  ^x  . 

(17) 
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Ai)  approximate  solution  to  the  contamiiution  problem  Lh  ol  'ained  b/ 
assuming  that  the  exposed  plane  passes  through  the  cloud  as  described 
at  some  fixed  time,  t.  The  aircraft  exposure  is  thus 

E(z,t)  a  /  X  dT  =  r  X  (x,z,t)  - —  = 

''o  b 

a 


Qo 


exp 


L 


(18) 


where  V_  ‘3  the  airspeed  of  the  exposed  plane:  or,  for  X  a  ut  ,  the 

tX 

position  of  the  center  of  gravity  of  the  cloud  from  its  initial  position. 


E(X,z)  -.j 


,1/ 


exp  !- 


^  2  '2-11 
^z 


(19) 


3.  VALUES  OF  PAH/UIETERS 

3. 1  Ground-level  Exposure,  Polri  or  Line  Source 

The  determination  of  suitable  values  for  u  ,  n,  C„  ,  and  C  for 

y  '  z 

this  part  of  the  problem  is  primarily  governed  by  the  fact  that  the  re¬ 
lease  is  at  a  relatively  high  level  in  the  atmosphere  whereas  estimates 
of  ground-level  exposures  are  required.  Ail  four  of  these  parameters 
vary  with  height  with  the  greatest  percentage  variation  probably  being 
that  of  the  wind  speed,  if .  The  parameters  also  have  horizontal  and 
time  variations,  but  these  variations  are  considered  secondary  in 
importance  to  the  vertical  variations  as  far  as  estimating  exposures  at 
the  ground  are  concerned. 

In  general  the  parameters  have  been  given  values  that  are 
averaged  over  the  depth  of  the  atmosphere  from  the  surface  to  the  re¬ 
lease  altitude,  h.  Several  heights  are  chosen  between  500  and  35, 000  ft 
for  this  purpose. 
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Values  of  u  were  derived  from  tabiilatlons  presented  by 

7 

Ciaitcher  ;  these  tabulations  provide  seasonal  means  of  wind  speed  and 
direction  as  well  as  time  varabillty  vs,  height  for  many  weather- 
observing  stations  In  the  Northern  Hemisphere.  Three  stations  are 
chosen  as  examples  to  Indicate  roughly  the  degree  of  variability  one 
might  expect  due  to  such  factors  as  geographical  location  and  season 
of  the  year.  The  wind  data  for  these  three  stations,  (Seoul,  Korea; 

San  Diego,  California;  and  San  Juan,  Puerto  Rico)  are  given  In  Table  1. 

The  directions  of  the  mean  winds  are  also  Important  and  should 
be  determined  for  the  region  of  interest  in  order  to  assess  downwind 
effects. 

For  planning  purposes,  estimates  of  extremes  of  exposure  ex¬ 
pected  at  the  ground  due  to  flight -Una  orientation  can  be  obtained  by 
using  values  of  9  a  go«  (Une  source  crosswind)  and  9  a  o«  (line 
source  parallel  to  mean  wind)  in  the  generalized  model. 

It  must  be  emphasized  that  the  actual  wind  observed  at  any  time 
will  generally  be  quite  different  from  the  mean  values  given  in  Table  1. 
For  a  given  mean  wind  direction  at  a  given  level  in  the  atmosphere, 
the  observed  wind  speed  can  be  expected  to  be  between  zero  and  the 
mean  speed  in  that  direction  25  percent  of  the  time. 

Additionally,  it  must  be  recognized  that  the  existence  of  wind 
shear  with  height  has  a  profound  effect  upon  exposure  estimates  at 
ground  level  in  that  shear  enhances  the  diffusing  power  of  the  atmos¬ 
phere.  This  meaits  that  the  model  tends  to  overestimate  the  exposure 
In  regions  of  appreciable  shear. 

Selection  of  suitable  values  of  n,  C„  ,  and  C,  must  take  Into 

y  z 

account  the  variation  of  these  parameters  with  thermal  stability  as 
well  as  with  height.  In  general,  n  increases  with  Increasing  thermal 
stability.  However,  available  observational  data  indicate  no  systematic 
variation  of  Cy  and  with  thermal  stability. 

Under  the  assumption  that  atmospheric  anisotropy  -'.fied  not  be 
accounted  for  by  assigning  values  to  an  Oy  and  an  ,  results  quoted 
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TABLE  1  (Cunt.) 

Saa  Diego,  Caliioriila  (32°  44'N;  117°  lO'W) 


HEIGHT  SEASON 

(thsds 

of  feet) _ 0-0.5 _ 5 

0-0.5  Winter  2  2 

Spring  3  3 

Summer  2  2 

Fall  2  2 

5  Winter  2 

Spring  3 

Summer  2 

FaU  2 

10  Winter 

Spring 
Summer 
Fall 

15  Winter 

Spring 
Summer 
Fall 

20  Winter 

Spring 
Summer 
Fall 

25  Winter 

Spring 
Summer 
Fall 

30  Winter 

Spring 
Summer 
Fail 

35  Winter 

Spring 
Summer 
Fall 


HEIGHT  (Thsds  of  feet) 


10 

15 

20 

25 

30 

35 

4 

6 

8 

10 

12 

14 

4 

5 

7 

9 

11 

12 

3 

3 

5 

5 

8 

7 

3 

3 

5 

8 

7 

8 

4 

6 

8 

11 

12 

14 

4 

5 

7 

9 

11 

12 

3 

4 

5 

5 

7 

8 

3 

3 

5 

6 

8 

9 

6 

9 

11 

13 

14 

15 

5 

7 

8 

11 

12 

13 

4 

5 

6 

7 

8 

8 

3 

4 

6 

7 

8 

9 

11 

13 

13 

17 

19 

8 

10 

12 

11 

15 

5 

0 

7 

8 

9 

5 

7 

8 

9 

10 

15 

17 

19 

20 

11 

14 

15 

18 

7 

8 

9 

10 

8 

10 

11 

12 

19 

21 

22 

18 

18 

19 

9 

10 

11 

11 

12 

13 

22 

24 

19 

20 

11 

12 

13 

14 

25 

21 

13 

15 

14 


TABLE  1  (Cont.) 


San  Juan,  Puerto  IB.co  (18°  28'N;  66°  Q7*W) 


HEIGHT 
(thsds 
of  feet) 

SEASON 

0-0.5 

5 

HEIGHT  (Thsds  of  feet) 

10  15  20  25  30 

35 

0-0.5 

Winter 

7 

7 

6 

5 

6 

9 

11 

14 

Spring 

5 

5 

4 

5 

6 

8 

10 

13 

Summer 

a 

8 

7 

7 

4 

5 

5 

6 

Fall 

5 

5 

4 

4 

4 

4 

4 

5 

5 

Winter 

7 

6 

5 

6 

9 

12 

14 

Spring 

5 

4 

5 

d 

8 

10 

13 

Summer 

3 

7 

7 

6 

5 

5 

6 

FaU 

5 

5 

4 

4 

4 

4 

5 

10 

Winter 

4 

4 

5 

7 

10 

12 

Spring 

2 

3 

4 

6 

8 

11 

Summer 

6 

6 

5 

4 

4 

5 

FaU 

4 

4 

3 

3 

3 

4 

15 

’.VL^tar 

3 

4 

7 

9 

11 

i 

Sprlro? 

4 

5 

7 

9 

12 

1 

Summer 

5 

4 

4 

4 

4 

FaU 

3 

3 

3 

3 

4 

20 

Winter 

5 

8 

10 

12 

Spring 

6 

8 

10 

13 

Summer 

3 

3 

3 

3 

FaU 

2 

2 

2 

3 

25 

Winter 

10 

13 

15 

SprLng 

10 

12 

15 

Summer 

2 

2 

3 

FaU 

2 

2 

3 

30 

Winter 

13 

17 

Spring 

14 

17 

Summer 

2 

3 

FaU 

2 

3 

35 

Winter 

19 

Spring 

19 

Summer 

3 

FaU 

1 

15 
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In  the  LSAEC  WASH  740  report,  by  Darad  and  Haugen,  “  and  by 

C  X  3 

Coi:valr“  are  essentially  consistent  with  Sutton's  recommendations 

for  the  variation  of  n  with  thermal  stability.  At  levels  In  the  atmos¬ 
phere  well  above  the  gradient  level,  however,  variations  In  thermal 
stability  are  relatively  small  and  a  value  of  n  slightly  gresder  than 
the  neutral  value  Is  probably  appropriate. 

Computations  of  from  data  published  by  WUkins^^  for  large- 
scale  diffusion,  assuming  an  n  value  of  0.25,  led  to  values  for 
that  show  a  slight  tendency  to  decrease  with  height  but  that  are  con¬ 
sistent  with  values  of  arising  from  several  other  sources;  foi 
example,  references  2,  9,  and  14. 

To  determine  suitable  values  of  C_  ,  three  recently  conducted 

m 

experiments  over  a  range  of  about  100  miles  were  analyzed  with  a  rate 
of  deposition  of  2  cm/sec  for  the  fluorescent  pigment  used  and  a  value 
of  n  of  0. 25.  The  release  altitude  for  these  experiments  was  about 
2000  ft.  Thermal  stability  ranged  from  slightly  stable  to  slightly  un¬ 
stable.  The  resulting  values  of  C,  are  essentially  consistent  with 
2  *  13 

Prairie  Grass  data  and  Sutton's  suggested  empirical  equation  for 

the  variation  of  "C"  wUh  height,  but  they  are  an  order  of  magnitude 
less  than  the  values  recommended  In  the  USAEC  WASH  740  report. 
However,  since  the  three  experiments  analyzed  are  directly  pertinent 
to  the  problem  of  line  sources.  It  Is  felt  that  the  analysis  of  these  data 
must  guide  the  choice  of  values  for  C_  until  further  experimental 
evidence  Is  available. 

The  values  recommended  for  n,  ,  and  as  a  function  of 
height  and  thermal  stability  are  summarized  In  TaUe  2.  Probable 
limits  of  these  values  are  also  tabulated  for  the  purpose  of  assessing 
the  degree  of  variability  one  nilght  expect  In  these  parameters.  It  Is 
emphasized  that  tMs  tabulation  of  values  Is  for  estimating  exposure  at 
the  ground  arising  from  release  of  material  at  a  given  altitude  h; 
that  Is,  the  values  represent  rough  averages  over  the  heights  Involved. 

For  the  small  particle  sizes  considered  In  this  study,  the  deposi¬ 
tion  of  material  from  a  cloud  to  the  ground  during  non- precipitation 
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x.iiiLE  2.  Values  of  a,  Cy,  and  vs.  heJgItt  of  release  aiul  thermal 
staoliity  to  use  In  models  for  estimating  exposure  at  tlie 
ground.  Uangea  of  the  values  are  iitdlcafed  in  parentheses. 


HEIGHT 

(ft) 

STABILITY 

a 

Lapse  0 

.40(0.2-0.8) 

0.05(0.01-0.1) 

II  »• 

0.15(0.0-0.  20) 

Neutral 

*•  tf 

0.  25(0.  20-0, 35) 

500 

Moderate  Inversion 

It  It 

tl  tl 

0.33(0.30-0.45) 

Strong  Inversion 

It  »l 

tt  tl 

0.50(0.45-0.80) 

Lapse 

II  II 

II  tt 

0.20(0.0  -0.20) 

Neutral 

It  It 

It  M 

0.  25  0.  20-0.  35 

5,000 

Moderate  Inversion 

It  tl 

fl  »l 

0.33(0.30-0.45 

Strong  inversion 

•  1  «• 

fl  It 

0.50(0.45-0.80) 

Lapse 

II  tl 

II  It 

0.20(0.0  -0.20) 

Neutral 

tl  II 

tl  II 

0.25(0.20-0. 35) 

10. 000 

Mo^leratft  Inversion 

II  II 

tt  It 

0.  33(0. 30-0.45) 

Strong  Inversion 

II  •• 

0.50(0.45-0. 80) 

15,000 

0.30(0.2-0.5) 

0.03(0.01-0.08)  0.30(0.20-0.35) 

20,000 

II  It 

tl  19 

II  tf 

25,000 

It  *1 

M  I' 

It  If 

30,000 

It  tt 

II  »» 

II  It 

35,000 

t#  tt 

II  It 

II  II 

conditions  Is  believed  to  be  limited  by  the  m';;leorologlcal  diffusion 
parameters  and  the  processes  of  Impaction  and  sticking  to  objects  on 
the  ground  rather  than  the  settling  of  particles  through  the  influence  of 
gravity.  Particles  brought  through  the  boundary  layer  to  the  ground  b>' 
the  turbulent  diffusion  process  deposit  on  vegetation  or  other  ground 
objects  by  Inertial  Impaction  and  diffusion  and  stick  by  electrostatic 
forces,  chemical  attraction,  or  other  means.  Investigations  of  this 
process  during  neutral  atmospheric  stability,  assuming  that  the  rate  of 
transfer  of  material  across  the  boundary  layer  is  equivalent  to  the  rate 
of  transfer  of  momentum  across  tho  boundary  layer.  Indicate  that  the 
velocity  of  deposition  should  vary  directly  v/itli  wind  speed.  There  Is 
also  reason  to  believe  that  the  deposition  velocity  should  change  with 
stability  In  the  lower  layers  of  the  atmosphere,  being  greater  durlrg: 
unstable  conditions  than  during  Inversion  conditions. 
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:it  IlniTford  i:ici  from 


Data  on  the  velccity  of  deposition  of 
Drltlsh  experiments  Indicate  a  value  of  2.7  cir./sec  is  appropriate. 
Calculations  of  the  deposition  velocity,  assui;  ing  that  the  flux  of  matter 
is  equal  to  the  flux  of  moinentuni  during  neutral  coiiditlons,  result  In 
about  the  same  value.  Measurement  of  the  velocity  of  depcaltion  of 
fission  products  resulting  from  arc  burning  of  uranium  has  indicated 
a  lower  deposition  velocity  of  about  0. 1  cm/sec.  The  difference  be¬ 
tween,  the  and  the  fission  product  aeroso’  values  is  believed  to 

be  due  primarily  to  a  net  difference  in  impaction,  retention,  or  adsorp¬ 
tion  efficiency  on  the  surface. 

For  the  model  proposed,  the  velocliy  of  deposition  la  assumed  to 
vai7  directly  with  the  wind  speed  near  the  ground  (say  2  meters  height) 
with  an  arbitrary  but  theoretically  tractable  allowance  made  for  the 
change  due  to  atmospheric  stability.  They  are  given  as  a  ratio  to  the 
wind  speed.  Separate  values  are  given  for  particulates  and^alogens 
based  on  the  experience  at  Hanford  and  English  experiments  for  the 
halogens  and  English  experiments  for  the  fine  particulates.  The 
noble  gases  do  not  appear  to  deposit  appreciably. 

In  an  analysis  of  rain-  scavenging  of  radioactive  particulate  matter 

Q 

from  the  atmosphere,  Greenfield  found  that  direct  interaction  of  rain¬ 
drops  and  particles  does  not  account  for  the  efficient  removal  of  material 
whose  diameter  Is  below  approximately  one  micron  diameter.  However, 
he  was  able  to  explain  the  removal  of  these  smaller  particles  by  allowing 
the  particles  to  mix  \vlth  the  water  cloud  before  the  rain  starts  so  that 
the  small  particles  that  are  scavenged  by  coagulation  are  then  placed  in 
a  position  to  be  more  efficiently  removed  l)y  the  rain. 

Our  present  knowledge  is  still  too  Inadequate  for  anytnlng  beyond 
a  suggestion  of  the  scavenging  processes  in  the  atmosphere  for  such 
finely  divided  material  as  envisaged  from  the  reactor.  Evidence  from 

world-wide  fallout  studies  suggest  that  the  very  small  fission  products 

:| 

become  attached  to  the  lutural  aerosol  particlr's  aiui  then  have  a  history 
in  th»-  atmosphere  slniilar  to  that  of  tne  host.  The  time  required  for 
•^uch  a  co'\lltlon  to  near  completion  Is  not  known,  but  one  would  not 
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TABLE  3.  Deposition  parameters  for  particulates  and  halogens  vs. 
thermal  staolUty.  After  Healy  (‘j;  p.  6) 


THERMAL  STADIUTi' 

Vg/U2 

Particulates* 

Halogens 

Strong  Inversion 

1.5  X  10-4 

2.4x10-3 

Moderate  Inversion 

2.2x  10-4 

3.4x10-3 

Neutral 

3  X  10*'^ 

4.6  X  10-3 

Unstable 

6  X  10-4 

8x10-3 

*  For  fission  products  attached  to  natural  aerosols  a  figure  of  about 
ona-thlrd  that  for  halogens  Is  suggested  rather  than  the  above 
figures  for  particulates.  This  figure  arises  from  qualitative  ap¬ 
praisal  of  world- wide  f^out  and  cannot  be  defended  qua.itltattvely 
at  this  time. 

expect  It  to  be  as  rapid  as  the  coagulation  with  larger  water  droplets 
because  of  the  relative  differences  In  mean  free  paths,  Junge^^  sug¬ 
gests  that  the  predominant  cause  of  the  modification  In  sir?  -distribution 
of  the  stratospheric  aerosols  on  the  way  down  through  the  troposphere 
is  the  repeated  cycle  of  condensation  and  evaporation  of  clouds,  a 
process  requiring  considerable  time. 

For  the  present  calculation  we  assumed  that  the  predominant 
mechanism  for  removing  these  small  particles  from  the  air  by  rain 
must  be  cloud-droplet- scavenging  coupled  with  .'ater  scavenging  of  the 
cloud  droplets  liy  larger  raindrops.  Effective  scavenglig  constants, 
that  Is,  the  time  required  for  the  number  of  particles  to  decrease  to 
1/e  of  the  Initial  value,  for  the  contaminant  activity  In  the  presence  of 
cloud  droplets  characteristic  of  stratus  and  cumulus  cloud  conditions 
were  computed  from  Greenfield's  curves  assuming  that  the  activity 
was  proportional  to  the  volume  of  the  particle.  The  cloud  character- 
l.dlcs  used  t)y  Greenfield  compered  favorably  with  later  data  reported 
by  aulm  Kampe  ard  ’.Velckmann'  and  were  not  altered.  Tlie  calculated 
times  required  for  one-half  the  activity  from  the  air  to  enter  the  cloud 
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TABLE  -1.  Parameters  for  cJoud -droplet  scaver:flr4j  of  particulate 
cloud  and  precipitation  scaventtAng  of  cloud  droplets. 


Height  of 
Release  (ft) 

Cloud  Type 

Rainfall  Kate 
( mm/hr) 

Scavenging 
Constant,  Jr 
(sec 'I) 

Elimination 

Constant, 

(sec'l) 

300  and 

5000 

Stratus 

0.5 

5  X  10'® 

2  X  lO'"' 

All  heights 
5,000-35,000 

Cumulus 

3.5 

a  X  10'® 

1  X  10'® 

droplets,  that  la,  the  >  values  called  for  in  the  model  after  mixing  of 
the  cloud  and  contaminant  particles,  are  given  In  Table  4. 

The  amount  of  activity  that  will  fall  out  as  rain  will  depend  upon 
the  time  of  mixing  of  the  contaminant  and  cloud  elements  and  the  rate 
at  which  the  cloud  droplets  are  swept  from  the  clnid  by  larger  rain- 
drops.  This  latter  aspect  of  the  problem  was  studied  by  Chamberlain'' 
uali:g  Langmuir's  theory  of  the  formation  of  raindrops  by  coalescen  .e 

4 

with  smaller  raindrops  and  Best's  relationship  between  raindrop  star 
and  rainfall  rate.  Results  of  Chamberlain's  study  applicable  to  this 
problem  are  also  summarized  In  Table  4. 

The  scavenging  of  fission  products  by  cloud  droplets  can  only 
occur  at  heights  at  wiiich  the  clouds  form.  The  heights  of  the  various 
types  of  clouds  vary  v/Uhln  wide  limits.  However,  scavenging  by 
stratus  clouds  should  be  limited  to  5000  ft  ard  below,  whereas  scaveng¬ 
ing  by  cumuliform  clouds  can  occur  at  any  height  between  2000  ft  and 
the  tropopause, 

3. 2  Exposure  to  Private  Aircraft,  Line  or  Point  Source 

Tor  this  part  of  the  problem,  only  the  determination  of  values  of 
u  follows  the  same  procedure  as  outlined  in  the  preceding  discussion. 
Seasonal  mean  values  for  the  three  stations  chosen  as  examples  of  the 
types  of  wind  distributions  to  expect  are  given  In  Table  1  (page  13). 

The  choice  of  values  for  C  ,  C  ,  and  n,  however,  requires  a 

y  * 
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somev.’hri.t  differont  procedure.  The  dJifuslon  problem  we  are  consid¬ 
ering  here  Is  generally  restricted  to  a  layer  of  the  atmosphere  that  Is 
not  particularly  deep  nor  close  to  the  surface.  If  the  layer  of  concern 
Is  well  above  the  gradient  wind  level  In  the  atmosphere,  that  is,  about 
10, 000  ft  or  more  above  the  surface,  one  wcnild  expect  the  turbulence 
to  be  nearly  Isotropic.  If  the  atmosphere  were  Incompressible  or  the 
mean  tropospheric  lapse  rate  neutral  rather  than  slightly  stable,  one 
could  use  theoretical  arguments  to  conclude  an  isotropic  field  of  turbu¬ 
lence  in  the  free  atmosphere.  In  any  event  It  is  difficult  to  conceive 
that  the  diffusion  coefficients,  and  C_,  can  differ  by  roughly  an 
order  of  magnitude  so  long  as  diffusion  through  the  atmospheric  friction 
layer  is  not  involved.  There  are  no  data  concerning  vertical  diffusion 
through  layers  of  the  free  atmosphere  that  can  be  brought  to  bear  on 
this  problem,  so  the  best  one  can  do  here  is  to  estimate  a  reasonable 

value  for  C_,  In  line  with  *he  discussion  outlined  above,  one  would 
z 

expect  C  to  be  slightly  less  than  C  in  the  mean,  although  at  times 
2  y 

it  could  conceivably  be  equal  to  or  greater  than  C^. 

Thus,  if  both  planes*  are  flying  at  10, 000  ft  or  more  above  the 
surface,  we  recommend  the  following  values  for  the  parameters;  , 
0,3;  C„,  0,1;  n,  0,30,  Variations  of  these  parameters  are  probably 
within  the  ranges:  0, 2  ''  Cy  <  0,5;  0,05  i  ^  0,2; 

0.20  C  n  ^  0,35. 

On  the  other  hand,  if  the  flight  altitude  of  either  plane  is  nearer 

5000  than  10, 000  ft  or  below  5000  ft,  the  values  of  C„,  C_  ,  and  n 

y  * 

given  in  Table  2  should  be  used. 

•I.  DISCUSSION 

Most  of  the  remarks  of  this  section  are  intended  to  provide  per¬ 
spective  for  use  in  evaluating  the  significance  of  exposure  levels 
estimated  from  the  foregoing  models.  In  general  the  approach  to  the 
problem  is  conventional  'with  methods  directly  obtained  from  the  Liter¬ 
ature  applied  where  possible.  It  is  hoped  that  individuals  using  this 
document  will  realize  the  uncertainties  and  the  lack  of  precision  of 
•A  privately-owned  or  commercial  and  a  nuclear- powered  aircraft) 
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many  cf  the  calculatioioi  techniques  used  nral  will  appreciate  the  roie  of 
the  "educated  guess,  "  both  in  deriving  the  modei  parameters  and  in 
Irgerpreting  the  results. 

Even  to  this  day  there  is  no  real  staiidard  set  of  raimbers  to  use 
in  the  diffusion  model  and  each  "expert"  will  use  his  own.  In  fact  none 
of  the  many  attempts  to  derive  expressions  for  the  diffusion  of  gases 
and  small  particulate  materials  in  the  atmosphere  can  be  rigorously 
defended  on  theoretical  grounds,  but  most  can  be  used  with  acceptable 
accuracy,  If  appropriate  values  of  the  diffusion  parameters  are  intro¬ 
duced  and  the  limitations  of  the  models  properly  assessed.  Corsequent- 
ly  there  Is  a  common  reluctance  among  atmospheric  physicists  to 
attach  quantitative  estimates  to  phenomena  so  vague  <*nd  uncertain  as 
depcultlon  and  ralnout  or  cloud  scavenging,  particularly  since  assign¬ 
ment  of  numerical  estimations  conveys  an  erroneous  Impression  of  the 
coiifidence  or  firmness  of  the  knowledge  constituting  the  basis  of  the 
estimate. 

The  atmospheric  diffusion  model  chosen  utilizes  the  fact  that  the 
damage  from  radioactive  materLols  Is  the  result  of  the  total  lidegrated 
concentration  or  the  product  of  the  concentration  and  time  and  not  of 
the  concentration  Itself.  Thus  the  results  are  expressed  In  Integrated 
exposure  regardless  of  the  time  cf  passage  of  the  material.  The  ex¬ 
posure  limits  to  oe  applied  la  calculatlj^  the  permissible  release  rates 
can  then  oe  expressed  In  terms  of  the  Integrated  concentration  at  the 
site  of  Interest,  with  the  amount  of  any  Isotope  In  existence  at  this  time 
being  corrected  for  the  decay  time  since  release. 

Care  must  be  used  in  summing  the  components  of  the  exposure  at 
various  points  In  the  environs  for  leternUiUiig  permissible  release 
rates.  The  external  exposure  to  humans  at  ground  level  will  arise 
from  material  remaliUng  in  the  cloud  and  that  deixjsited  on  the  ground 
so  tliat  the  change  In  size  of  the  cloud  and  deposition  pattern  with  dis¬ 
tance  must  be  taken  Uito  account,  burther.  Isotopic  fi  actlonatlon  of 
both  the  material  In  tiie  cloud  and  that  depo  dtod  may  occur  at  dowinvlrd 
points  due  to  diilcrciU  deposition  velocities  of  the  Isotopes,  yVll  of  these 


factors  must  be  properly  accounted  for  in  the  model. 

The  diffusion  model  Is  only  expected  to  provide  an  average  be¬ 
havior  of  the  contaminants.  For  essentially  instantaneous  releases,  a 
large  variety  of  behaviors  Is  possible.  Discrete  puffs  are  Influenced 
completely  by  the  eddy  structure  In  which  they  are  embedded.  However 
this  structure  may  be  changing  constantly  so  that  the  variety  of  patterns 
Is  almost  Infinite. 

As  la  the  practice  wUh  calculations  of  this  nature,  the  dilution  or 
exposure  Is  calculated  for  a  specific  meteorological  condition.  However, 
meteorological  conditions  are  characterized  by  change  so  that  the  proba¬ 
bility  of  encountering  protracted  periods  of  a  specific  meteorological 
condition  la  rather  low,  especially  In  the  lower  atmosphere.  Thus  the 
"steady  state"  estimates  made  for  large  travel  times  may  become  In¬ 
creasingly  unrealistic.  The  mean  wind  speeds  given  in  Table  2  (page 
17)  are  averaged  over  the  space  governed  by  the  radiosonde  flight 
column  to  35,000  ft  from  a  given  station.  In  the  aircraft  flight  problem 
not  only  the  wliuls  in  the  area  of  flight  but  the  winds  in  the  area  affected 
by  the  contaminants  should  be  characterized  In  terms  of  space  and  time 
variability.  Such  a  climatic  summary  for  the  specific  areas  Is  required 
to  adequately  appraise  the  environmental  problem.  Tlien  one  can 
assess  the  probability  of  certain  levels  of  environmental  effects  based 
on  the  climatic  summaries.  Wind  velocity  shear  Ln  the  vertical  Is  an 
important  factor  in  dispersing  materials  In  the  atmosphere,  particularly 
over  long  travel  times.  In  the  lower  5COO  ft  of  the  atmosphere  the 
average  veering  is  20**  to  90”.  To  a  first  approximation,  one  may  as¬ 
sume  that  this  veering  Is  evenly  distributed  throughout  the  layer. 

Another  factor  not  specifically  Irxluded  in  the  model  Is  the  effect 
of  large-scale  vertical  motions  In  the  atmosphere.  However  It  Is  felt 
that  the  effect  of  this  factor  Is  sufficiently  allowed  for  the  rar,ge  of 
C_  appearing  in  Table  2. 

Altitudes  as  boundaries  for  certain  physical  processes  In  the  at¬ 
mosphere  are  not  liUended  to  be  rigidly  defined.  In  deriving  the  values 
of  parameters  and  altitudes  of  appllcablllly,  some  space  and  time 
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average  thought  to  be  reasonable  was  used  to  provide  an  internally  con¬ 
sistent  model.  In  some  cases,  only  an  educated  guess  was  available. 
Particularly  dilficult  or  Impossible  to  incorporate  in  a  general  model 
is  the  effect  of  changes  of  thermal  stability  with  altitude.  Vertical 
temperature  soundinga  commonly  show  layers  of  stable  and  unstable 
air  at  various  altitudes.  Vertical  diffusion  of  contaminants  released 
between  stable  layers  may  be  limited  almost  as  If  bound  by  horizontal 
barilers  above  and  below.  An  analysis  of  temperature  soundings  In 
the  region  of  Interest  would  Indicate  the  severity  of  such  a  problem. 

In  coastal  regions  a  land  and  sea  breeze  may  have  to  be  accounted  for 
In  the  analysis.  The  vertical  extent  of  sea  breezes  is  normally  about 
3000  ft.  the  land  breeze  about  1500  ft,  and  the  overall  land- sea  breeze 
system  In  the  vertical  extends  to  about  10. 000  ft. 

Noteworthy  among  the  meteorological  situations  not  adequately 
covered  by  a  general  model  Is  one  such  as  the  capping  inversion  In  the 
Los  Angeles  area.  Dispersion  of  pollutants  released  In  the  layer  be¬ 
neath  the  Inversion  Is  largely  confined  to  this  layer,  while  horizontal 
dispersion  is  largely  confined  to  the  areal  extent  of  the  land- sea  breeze 
system.  A  climatic  analysis  of  the  region  would  Indicate  the  frequency 
of  such  adverse  meteorological  conditions. 

Aiuther  example  of  such  a  region  Is  found  near  the  trade-wind 
belt  In  tropical  latitudes  (roughly  5’N  to  25^N  In  the  Northern 
Hemisphere).  These  so-called  trade-wind  Inversions  are  particularly 
persistent  over  ocean  surfaces  or  small  Islands. 

If,  on  the  other  hand,  releases  are  consistently  made  well  above 
the  capping  Inversion  In  regions  such  as  these.  It  Is  to  be  expected  that 
exposure  levels  at  the  ground  would  be  reduced  from  those  estimated 
using  the  models. 

Another  type  uf  meteorological  situation  somewhat  similar  to 
that  just  discussed  from  the  poUd  of  view  of  estimating  exposure  levels 
at  the  ground  Is  the  so-called  polar  outbreak.  In  this  case  a  mass  of 
cold  atr  travelir4{  behind  a  cold  front  Is  normally  capped  by  a  strong 
Inversion,  Once  again,  roleases  beneath  the  cap  of  the  polar  air  mass 
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would  Increase  the  hazard  problem  at  the  surface  for  relatively  large 
areas  siiu:e  a  strong  polar  outbreak  characteristically  moves  over  hun¬ 
dreds  of  miles  without  significant  change. 

It  should  also  be  noted  that  the  model  does  not  account  for  changes 
in  the  stability  of  the  air  in  the  lowest  5000  ft  of  the  atmosphere  as  the 
air  travels  from  land  to  water  or  vice  versa. 

The  atmospheric  diffusion  of  the  fission  products  that  come  in 
contact  with  the  ground,  vegetation,  buUdlngs,  etc. ,  is  complicated  by 
the  loss  due  to  deposition,  which  tends  to  deplete  the  layer  of  the  plume 
in  contact  with  the  ground,  leaving  an  additional  radioactive  exposure 
from  the  nuterials  so  deposited.  The  key  parameter  In  the  deposition 
problem  is  the  deposition  coefficient  for  which  reliable  values  must 
come  from  experimental  measurements.  Even  so,  the  deposition 
phenomenon  is  believed  to  be  so  complicated  by  interactions  among 
diffusion  parameters,  surface  characteristics,  contaminant  interactions, 
carrier-host  problems,  retention  efficiencies,  and  chemical  nature  of 
contaminants  and  surfaces  that  interpretation  of  results  is  very  difficult. 

Cloud- dropli-t  scavenging  and  raindrop  scavenging  are  simplified 
in  the  model.  It  has  proved  suitable  in  meteorological  practice  to 
characterize  a  rainfall  by  rainfall  intensity  and  the  average  diameter 
uf  drops  and  their  terminal  velocities  (Oest  ).  The  computed  values 
correspond  only  approximately  to  the  actual  conditions,  as  the  size  of 
the  drops  always  shows  a  spectral  distribution.  The  curve  of  Dest  con¬ 
tains  mean  values  which,  in  single  cases,  show  remarkable  differences, 
with  the  extreme  range  in  rainfall  rate  for  a  given  mean  drop  size  about 
a  factor  of  two. 
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Sep  1953.  (SECRET  Report) 

.No.  48.  Oparatioa  t'pabot>Knotbola  Project  1.3.  Free  Air  Atomic  Blast  Preaaara  Measuremaats.  Re« 
vised  Report  (V),  ,V.  A.  Haskell  and  R.  If.  Brubaker,  .Sou  fOSJ.  (SECRET/ RESTRICTED 
DATA  Report) 

No.  49.  Maiimom  Humidity  ia  Eapinceriop  Desifm.  -A.  Sissenuine,  Oel  1953, 

No.  so.  ProbeUa  lea  NlaaJ  l.ocatioas  in  tbe  Arriie  Basin,  Janaarv  19.34,  .4.  P.  Crary  and  I,  Browne, 
May  1954. 

No.  51.  Inveaii;;atioa  of  TR.AC  for  .Active  Air  De^enaa  Purposes  (L),  G.  9.  foree,  R.  Penndorf,  F.  G, 
Plank  and  B.  H.  Grossman,  Dee  1953.  (SECRET  RESTRICTED  DATA  Report) 

No.  52.  Radio  Noise  Emissions  During  Tbermoauclcar  Reactiona  (VK  T,  J.  Keneekett,  Jun  1954. 
(CO.SFlDt.STIAL  Report) 

No.  .33.  A  Metbod  of  Correetiatt  Tabulated  Rawinaonda  Rind  Speeds  (or  Curvature  of  the  Kartb.  R, 
Levitan,  Jun  19 jt. 

No.  54.  .A  Proposed  Radar  Storm  Wamioa  .Service  For  Army  Combat  Operations.  If.  C.  //.  LteJa,  Aug 

1954. 

'lo,  3.3.  .A  Comparison  of  Altitude  Corrections  for  Blast  Overpressure  (V),  .V.  .1.  Haskell,  Sep  1954, 
(SECRET  Report) 

No.  36.  Attenuatintt  Effects  of  Atmospheric  i,i‘]uiJ  Rater  on  Peak  Ovrrprearurrt  from  Blast  Waves  (V), 
H,  P.  Cauvin,  ].  H.  Healy  and  V.  .4.  Bennett,  Ort  I95t.  (SECRET  Report) 
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\o.  '7.  Sin  Ispccd  Profile,  Sinilshedf,  and  (jiisis  for  OeniKn  cf  (Guidance  Syalema  for  Vertical  Hiainf; 
\ir  'i -iiii  let,  \.  N  iaenuioe,  .Vot/  J9St. 

No.  .mI.  Phe  Snppreasion  of  Aircraft  Fxhaust  Praila,  h.  AnJerton,  Sov  1954. 

No.  '>9.  Preliminarv  Keport  on  lha  Attenuation  of  Thetmal  Hadiatinn  Front  Atomic  or  Pfiennonuclear 

iVeapona  tUI,  H.  tf.  Chapman  and  W.  //.  Seavev^Mov  1954.  (SECKET/RCSTRICTED  DATA  Re^ 
part) 

No.  60.  ilaigbt  turora  in  a  Haario  Syatem,  R.  Laviton,  Dee  1954^ 

No.  61.  Metaaroionical  Aapacta  of  Conatanl  Laval  Balloon  OpatMioaa  fL'J^  f.  K.  ffidger,  Jr.  tt  at,  D*e 

1954.  (SECRET  Report) 

No.  62.  V'ariaiioaa  in  (^aomatric  Haight  of  30  to  60  Tlioiiaaad  .‘‘‘oot  Preaanra-Altitudaa  (UX  Af. 
Siteenwine,  A,  E.  Cole  and  9.  Baginaky,  Dee  1954.  (COSFIDENTIAL  Report) 

No.  63.  lirview  of  Tima  and  Spara  Kind  Kluctualiona  Applicabfa  to  Convantional  Balliaiie  0e)ar> 
ininaiinna,  V,  Baginsky,  \.  Sittenwine,  B.  Davidton  and  H.  Lettau,  Dee  1954, 

No.  6t.  (^loudinaaa  Abova  20,000  Feat  (nr  Cartnin  Stallar  Naviuntioo  Problana  (OX  A.  E.  Cole,  Jaa 

1955,  (SECRET  Report) 

So.  65,  Tha  Kaaaibiliiv  of  tha  Idaniificatioa  of  Hail  and  Savara  Sionaa,  D.  Atlai  and  R.  DonaJdaon, 

Jan  1955, 

No,  66.  Rala  of  Rainfall  Frei{aaaeiea  Ovar  Salarlad  Air  Routaa  aad  Oaaiiaationa  (UX  A,  E.  Cole  and 
\.  Sieeemeine,  dar  1955.  (SECRET  Report) 

No.  67.  Stima  Coaaidarationa  oa  tha  Modaliofi  of  Craiarinfi  Pbaaoaiaaa  in  Fartbff/A  Af.  A.  fjaekell,  Apr 
1955.  (SECRET.  RESTRICTED  DATA  Report) 

So.  6A  Tha  Preparation  of  Kvtanded  Foracaaia  of  tha  Preaaara  Heiaht  Oiatribnlion  in  tha  Frea  Atmon* 
phere  Over  Noitli  \mertra  bv  I'aa  of  Kmpirical  Indueaca  Functiona.  fi.  ’J,  Shite,  Way  1955. 

No.  69.  C.ild  Keathar  Fffeci  on  iM2  Launchinn  Pemoanal  (U),  \.  Sittenuine,  Jun  1955,  (SECRET  At* 
port) 

No.  70,  Vtmrapheric  Preaaura  Pulaa  Mraaoremenia,  Operation  Caalla  (L'X  E.  A.  Flauraud,  Aug  1955. 
(SECRET  RESTRICTED  DATA  Report) 

Vo.  71.  Ilefraetioa  of  Shock  Ka»ea  in  ih*  Atmoaphere  (L'X  A,  f.  Ilaekelt,  Ang  1955  (SECRET  Report) 
Vo,  72.  Kind  •ariabililv  aa  .i  Fuaction  of  Time  at  Muroc,  California,  B.  Singer,  Sep  1955. 

No.  73.  The  Atmoaphere,  .V,  C.  ('•erton,  i  ‘p  1955. 

So.  71.  kreal  \  ariatioa  of  (.'eiiinn  Heiaht  (i),  S,  Bagineky  and  A.  E.  Cote,  Oet  1955.  (COSFIDESTJAL 
Report) 

So.  73.  An  Objective  .Svaiem  for  Preparing  Operational  Weather  Foracaata,  I.  A.  Lund  and  E.  f,  takl, 

.Vot'  1955. 

So.  76.  Ihe  Practical  Vapecia  of  Tropical  Metenraloay,  C.  E.  Palmer,  C.  If.  9ise,  L.  J,  Stempton  and 
C.  II,  Duncan,  Sep  1955. 

Vo.  77.  Ilemola  Deterniinalinn  of  Soil  Traffieability  by  Aerial  Penetrometer,  C,  Molineux,  Oet  1955. 

So.  7A.  Kffecta  of  the  Primary  Coamic  Radiaiion  oa  Matiar,  11.  0.  Curtis,  Jan  1956. 

Vo.  79.  rropoaphene  '  oriationn  of  Refracti v  .'ndex  at  Mict'wava  Freoaentiea,  C.  F,  Campen  and  .4.  E. 
t.alf,  Oct  1953. 

No.  SO.  V  Procram  to  Pent  Skill  in  Teiminal  ‘  .ireceaiinji,  /.  I.  Crtngorten,  I.  .4.  Lund  and  if.  A.  Viller, 

Jun  i  1 55. 

V.i.  111.  I'itir'-me  Vimoapherea  and  Hallialic  .  <.itir*,  .V.  .Sittenitine  and  A.  E.  Cnie,  Jul  1955. 

No.  82.  Hot  ifi.in.il  h re-pienciee  and  Vbsorption  Coefficient*  of  \:moapheric  Caeee,  S.  .V.  Chash  and 
It.  IK  t-.litar.ls.  Mar  l')56. 

No.  'll.  lon'Upheric  h.ffrrln  on  Poailionina  of  Vehicle*  ot  Mich  Altitude*,  if.  (' filter  and  T,  J. 

Rcncxhca,  Mar  1956. 

V.I.  .1,.  I're-I'rouah  Sinter  Precipit.iti.m  Forecaatina.  P.  S.  Faiae,  /’eS  (957, 
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No.  8S.  Cecmagnetic  Field  Extrapoiatioe  Teciiaiquea  —  Aa  Evaluation  of  the  Poitaon  Inte^rai  for  a 
Plana  (VK  /•  F.  UcClay  <ini  P.  Fougert.  Fei  1957.  (SECRET  ReporO 
No.  86.  'rha  ARDC  Model  Abcoaphcra,  1956,  R.  A.  .Vtuaaroad  7.  5.  RipUy,  Dee  1956. 

No.  87.  Aa  £atima:e  of  the  Maxiainni  R*>R*  '>•  Detectability  of  Sai  mic  Sii^aala,  A.  A.  Haektll,  Uar 
1957. 

No.  88.  Sooia  Coacapta  for  ProdictiaD  .Naelaar  Cratar  Siio  (U)^  F.  A.  Crowley,  Fek  1957.  (SECRET/ 
RESTRICTED  DATA  Reocet) 

No.  89.  Upper  ffiad  Rapreaeatatioo  aad  Flight  Plnaaiai^  /. Griegortie  Mi  r  1957. 

No.  90.  Rafleetioa  of  Point  Soarco  Radiatios  From  a  Lambort  Plaaa  Onto  a  Plaaa  Ree*>ear,  A.  f. 

Cueee,  Jel  1957. 

No.  91.  Tba  Va.iatioaa  of  Atmoapbarie  Traaaniaaieity  and  Cloud  Height  at  NawaHt,  T,  0.  Haig,  aed 
t.  C.  Morton,  III,  Jan  195S. 

No.  92.  Collaetioa  of  Aaromagaede  laforaatioa  For  Goidaaeo  aad  .Navigation  (Vk  R.  Hnukinton,  B. 
Skmmon,  R.  Brsreton  and  J.  McClay,  .fag  1957.  (SECP.ET  Report) 

,No.  93.  Tba  Aeenmcy  of  Wind  Deteraiiaation  From  tba  Track  of  a  Falling  Objaci,  F.  Lolly  and  R. 
LeoUon,  Mar  1959. 

No.  94.  Cadmadeg  Soil  .Moiatara  aad  Tmelionability  Coaditicna  for  Simtagie  Plasaing  (V)y  Part  1  • 
Gaaaml  niatkod,  and  Part  3  •  Applicaiiona  aad  iaiarpratadoaa,  C.  7.  Tkomtkwoiu,  /.  R. 

Matket,  D.  8.  Carter  and  C.  E.  IhlinenM,  Mer  1959  (Inelaeeified  Report).  Part  3  •  Avaiaga 
aoil  moiataro  aad  traedoaabilily  conditioaa  ia  Poland  (V),  D.  B.  Carter  and  C.  E.  MoHnenn,  dag 
1959  ((XJNFIDESTIAL  Reportji  Part  4*  Avorago  aoil  moiatara  aad  trncdoaabilitv  coaditioaa 
ia  Yngoalavia  (DJ,  D.  B.  Carter  and  C  E.  Molinene,  Ma.  1959  (CONFIDESTIAL  Report) 

No.  9S.  Rind  Spaeda  at  50,000  to  1000,000  Feet  and  a  Rrlatad  Balloon  Platform  Dtaign  Problem  (U),  .V. 
Dvoekin  and  S.  Sieeentnine,  jitl  1957.  (SECRET  Report) 

No.  96.  Davalopmaal  of  Miaaila  Oaaign  tiad  Prontaa  for  Patritk  AFB,  .V.  Sireemtnne,  Mar  1959. 

No.  97.  Cloud  Baaa  Dalaetioa  by  Airborne  Radar,  R.  J.  Donoldton,  Jr.,  Mar  1959. 

No.  98.  Maoa  Fraa  Air  Gravity  .Aaomaliea,  Gaoid  Coatoor  Cervao,  aad  tba  Avorago  Defleetiona  of  tba 

Vertical  (OX  r.  A.  HeUknnen,  U.  A.  VotiUendO.  V.  Billumt,  Mar  1959.  (COSFIDESTIAL  Rf 
port) 

No.  99.  E'-'alaalioa  of  A.N/CMO-3  Wind  Sbeor  Data  lor  Dovalepmoat  of  Miaaila  Daaiga  Criteria,  S. 

'  joekin  and  .V.  Staaaaaraa,  Apr  1959. 

■NoJOC.  .A  Pbenomenological  Tbaory  of  tba  Scaliag  of  Fireball  Miaimem  Radiant  latenaity  with  Yield  and 
Altitada  (UX  H,  K.  Sen,  Apr  1958.  (SECRET  Report) 

N.v.lOl.  Cvataation  of  Satellite  Oboarving  Network  for  Project  “Space  Track",  C.  R.  Mietaika  and  H.  0. 
Curtis,  Jun  1959. 

.No.  102.  Aa  Oparatioaal  Syatam  to  Mraauca,  Compote,  and  Praaaat  Approach  Viaibilily  lafomiatioa,  T.  0. 
Haig  and  W.  C.  Morton,  Hi,  Jun  1959. 

No.l0\  llararda  of  l.ighlaiag  Oiacbargo  to  .Aircraft.  G.  A.  Faueker  and  H,  0.  Cui:<t,  Aug  1959. 

No.104.  Conuoil  Prodiclion  and  Prevention  (I'X  C.  S.  Donnie,  C.  £.  Anderson,  S.  J.  Birstein  and  B.  A. 
Silverman,  Aug  1959.  (SECRET  Report) 

No. 105.  Meihoda  o.  Artificial  Fog  Oiaperaal  and  Their  Kvaluatioa,  C.  E.  Junge,  Sep  1959. 

No. 106.  Themal  Tachaiqnaa  for  Oiiaipating  Fog  From  .Aircraft  Runwaya,  C.  5.  Downie  and  R.  B.  Smith, 
Sep  1958. 

No. 107.  Accuracy  of  RDF'  Poaition  Fivca  in  Tricking  Cooaimt>l.evel  Ballooaa,  K.  C.  Giles  and  R.  E. 
Peterson,  edited  by  T.  K.  tiJger,  Jr.,  Oct  1958. 

No.lOfl.  The  Effect  of  ifind  Errora  on  SAGE-Guided  Intercepta  ff/Jl  E.  M.  Darling,  Jr.  at.d  C.  D.  V-tt, 

Oet  1959  (COSFIDESTIAL  Report) 

.No.  109  Behavior  of  Atmoaphenc  Denaiiv  Profilea,  .V.  Sirrcnu'inr,  7.  .S,  Ripley  and  .f.  E.  Cole,  Dec  1958. 
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No.!  10.  Mo^neiic  Ooterniinaiion  of  Spaca  Vebicle  Allitude ].  h.  \lc(,lay  uud  I‘,  h.  ''ou^ere,  ^la’ 

1959.  (SkCKKT  Report) 

No. III.  Final  Report  on  Exbaaal  Trail  Phyaica:  Project  76.10,  T.iek  76108  (  f).  If.  //.  .WcAennc,  unJ 

il.  0.  Curtia,  Jul  1959.  (SECKKT  Report) 

No.  112.  Accaracy  oi  \laaB  Moatlily  Gaoatrophic  Riad  Veciora  aa  a  Function  of  Station  NetwoHt  iJen- 
■•ity,  H.  A.  Saimtl(^  ]m  1959. 

Na.ll3.  Aa  Catiaiaia  of  iba  Stf«n((tli  of  dia  Ae««atie  Sigar  l  Caaamtcd  br  aa  ICIIM  .None  ^oaa  Reaatry 
,'f/A  N.  A.  H'uMl,  Ah  J959.  (CONFIDENTIAL  tiepanj 

No.llA.  Tfca  Rola  of  Radiatioa  ia  Shoefc  PropoKOtioa  widi  Applicatieaa  to  Altitado  oad  Yield  ScallaK  of 

Noclo«r  Firobidla  (U),  H.  K.  Son  mA  A.  f.  &'««aa.  Sep  1959.  (SKCUET/ULSTRIC1..U  DATA  Repent 

.No.lIS.  .AROC  .Modal  AtmoopAaro,  1959,  R,  A.  IlioMmer,  K.  S.  9.  Ckampioit  and  R.  /,.  Pond,  iun  1959, 

No.116.  llofiaoMoata  in  Utilitaiion  of  Contour  Chario  fot  Cliniatirully  Specified  Wind  Profilea,  A.  E. 
Cole,  Oct  1959. 

N".117.  Dtoiffa  Vlad  Profilao  From  Japaaaao  Relay  Souadiag  Dau.  N.  Sittenwine,  U.  T.  Mulkttn, 
and  h.  A,  SalmeUt,  Dee  1959. 

N0.U8.  Uiliuor  Applicatioaa  of  Supercooled  Cloud  aod  Fog  Diaotpniion,  C.  S.  Donme,  and  B.  A. 
SUvenmm,  Dee  1959. 

No.  119.  Fueler  Auelyeiu  aod  Siep«>iae  Regroeeioa  Applied  to  lite  2A.llour  IVedicliou  of  SOOhuIi  Vindo, 
Taoipamiurua.  and  lloighlu  Over  a  Siiaoi  Arau  (VK  E.  J.  Aeberi,  I.  A.  I.utd,  A.  Thomaeell,  Jr,. 
and  J.  /.  Patruokaa.  fek  I960.  (CONFIDENTIAL  HepeH) 

No.l2a  Aa  GatiMle  of  Prerjpiiable  Raler  Aioog  llieii.Altiiade  Hay  P.itha,  Ifumv  Cutniek,  Mar  196’.. 

No.i21.  .Vnairtiag  and  Foreeaaiiag  .\Oi|aorolopical  Coodiiioae  ia  tbe  Upper  IVopoapliero  and  (attver 
.ttratooplMro,  R.  M.  f.ndUeh  and  C.  3,  Me  Learn,  .Apr  I960. 

No.122.  Aaalyoia  and  IVedicliea  of  the  S00>mb  Surface  ia  a  Silent  Area,  (C),  E.  A,  Aekert,  May  I960 
(CONFIDENTIAL  Repertk 


